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Abstract 

To help understand the centre dominance picture of confinement, we look at Wilson loop distri- 
butions in pure SU(2) lattice gauge theory. A strong coupling approximation for the distribution 
is developed to use for comparisons. We perform a Fourier expansion of the distribution: centre 
dominance here corresponds to suppression of odd terms beyond the first. The Fourier terms 
correspond to SU(2) representations; hence Casimir scaling behaviour leads to centre dominance. 
We examine the positive plaquette model, where only thick vortices are present. We show that a 
simple picture of random, non-interacting centre vortices gives a string tension about 3/4 of the 
measured value. Finally, we attempt to limit confusion about the adjoint representation. 
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Figure 1: Monte Carlo distribution of Wilson loops on a 12 4 lattice at (3 — 2.5, showing also the 
limit of the distribution for large loops. The area under each curve is normalised to unity. On this 
scale errors are negligible. 

1 Introduction 

Recently the idea, originally proposed some time ago [jl], |[ [| , that confinement in gauge theories 
can be considered as an effect related to the centre Z(N) of the gauge group SU(iV) has been 
undergoing a revival 0, 0, 0, M. Not the least of this has been the demonstration H that if 
one replaces the values of Wilson loops with their signs alone (i.e. the centre Z(2)) the heavy 
quark potential in SU(2), which is now very accurately known can be essentially completely 
reproduced. This is a spectacular — and gauge invariant — demonstration that something is right 
about this hypothesis; it is, however, far from a demonstration that centre vortices are responsible 
for confinement. 

The key part of the argument for Z(2) to be an ingredient in confinement is that there are 
'thick' centre vortices, associated with the quotient group SU(2)/Z(2) = SO(3), which pierce 
Wilson loops; their physical effect depends on how many vortices, modulo 2, pierce a given area. 
There are also 'thin' vortices associated with Z(2), which appear as chains of negative plaquettes; 
to create these requires an action proportional to the number of flipped plaquettes, so these will 
not survive in the continuum limit. By contrast, the thick vortices show up only in larger Wilson 
loops, which may be negative while still surrounding only positive plaquettes, and do survive in 
the continuum limit; they are topological in nature, related to the fact that SU(2)/Z(2) is not 
simply connected and indeed has a Z(2) homotopy. (For simplicity, we have here ignored hybrid 
vortices, which combine the two effects.) We shall return to the distinction between Z(2) and 
SU(2)/Z(2) effects towards the end of the article because, as emphasised in ref. where this 
mechanism is described in more detail, it is important and has caused much confusion. 

We have specialised to SU(iV) for TV = 2 since the arguments are expected to extend to higher 
N via the centre Z(TV), though of course this needs to be checked explicitly. We are also ignoring 
the effects of fcrmions in the vacuum. 
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Because Z(2) commutes with all elements of the gauge group, the sign of any Wilson loop in 
the fundamental representation of SU(2) must be determined by the combined effect of thin and 
thick vortices: the total of these objects passing through the loop determines whether the loop 
is positive (even number of vortices) or negative (odd number). In an attempt to see how this 
appears, we can look at the distribution of the trace of the Wilson loop, as generated by a Monte 
Carlo simulation. Note we are looking at the distribution of individual, not average, Wilson loops. 
This was first looked at, in the case of the plaquette, in refs. |], [To| , where it is referred to as the 
'spectral density', and the point was made that it contains information about all representations. 
One such set, taken from 10,000 configurations of a 12 4 lattice at j3 = 2.5, is shown in fig. [IJ It is 
clear that it changes smoothly from the highly asymmetric shape of the plaquette action to some 
limiting distribution for large loops. It is less clear how the Wilson loop expectation value (Wo), 
which is the first moment of the distribution, behaves; nor is it clear where the special role of the 
sign of the loop comes from. 

This paper attempts to shed some light on these matters. As the behaviour of the distribution 
is largely unfamiliar territory, we first perform a calculation using a simple approximation: we are 
not able to calculate the raw distribution, but we are able to see how it develops for larger loops. 
The expectation value in this approximation has the leading order strong coupling behaviour. 

We shall proceed as follows: first (section 2) we shall give a few basic home truths about 
Wilson loops. Next (section 3), we shall look at the approximation that all loops are uncorrelated, 
corresponding to strong coupling; in section 4 we discuss the centre dominance picture in the 
same spirit, where we also explicitly show the connection between the shape of the distribution, 
irreducible representations of SU(2) and the phenomenon of centre dominance. Then we make some 
remarks about centre-projected vortices, the positive plaquette model, the adjoint representation, 
and show how a gas of non-interacting Z(2) vortices gives exact area law confinement. In the last 
section we summarise these results concisely. 

Section 3 can safely be missed out by anyone not interested in the details of the distributions of 
uncorrelated Wilson loops. Of the formalism in that section we shall only use the formula relating 
the distributions for loops of areas A and 2 A, eqn. (|l2|), and its extension to arbitrary multiples 
of A. The main physical discussion is in section 4. 

2 Wilson loop distributions 

We shall start by assuming the standard Wilson lattice gauge theory in SU(2). We parametrise 
an open Wilson loop (i.e. the SU(2) element representing the loop, where no trace has been taken) 
as 

W = W a l 2 + i<7.W, Wq + wl + w£ + w£ = 1. (1) 

We now consider the distribution, the word being used in the sense of a Monte Carlo calculation: 
we have a large number of configurations, and consider the spread of values for all Wilson loops 
of a given shape on all lattices. Of course, we could equally well talk in terms of contributions 
to the path integral, where, because of the importance sampling built into the Monte Carlo 
simulations, the exponential factor is included and there is a uniform measure; this is entirely 
equivalent. However, the statistical language fits well here where we use actual Monte Carlo data 
for comparison. The distribution contains all the gauge invariant information about the loops, 
including the expectation values of the loops in all representations. Indeed, in a quantum theory 
it is natural to consider this distribution as the basic quantity. 

The first matter of interest is the limiting large loop distribution apparent in fig. ^. This can 
easily be calculated. Given the short range behaviour of the force, for sufficiently large loops 
different parts of the loop are physically unconnected with one another. The lack of any overall 
correlation means the loops are random: the distribution corresponds to a random walk by Wi 
over the 3-sphere of the gauge manifold, so we simply need to calculate the fraction of the surface 
available at a given Wq with the other co-ordinates fixed. The probability that Wq lies in a certain 
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range is 



p l (W )dW =S ds /S, (2) 

where here, as throughout, we use the symbol p for the distribution and the suffix I indicates 
the limiting distribution, Sq is the 2-surface for a given Wq, while between Wq and Wq + dW a 
distance dso is covered on the surface of the 3-sphere with the other co-ordinates held constant, 
and S is the total surface area of the 3-sphere. Basic geometry gives 

nm)dWo . MLl*a_j^ _ 2 (1 _ <y n iWa , (3) 

This is the lowest curve in fig. [j] and clearly fits the part. Indeed, loops larger than 4x4 are so 
close to this that they have been omitted for clarity. Later, we will see that this curve corresponds 
to the identity representation of SU(2); it will then be obvious that convergence to this limit is 
inevitable provided that all j = 1/2 and higher representation Wilson loops decay monotonically. 

The loop distribution in figure [l| is smooth and we would expect to be able to expand it in 
a Fourier series. We will pick the units Wq = sin^, where — ir < <fi < it so that the terms are 
orthogonal on the range of values of Wq. However, our expansion is not quite standard; we use 
our freedom to expand the even and odd parts separately to write, 

Pa(W ) = pa (sin — ) = ^(a„cos(n- -)<j) + b n siring) (4) 

n— 1 ~ 

as the distribution of Wq at loop area A. The difference from the usual Fourier series is that this 
form satisfies the boundary condition that the function vanishes at cf> = ±ir. The basis corresponds 
to the eigenf unctions of the self-adjoint linear equation 



d 2 p 

d(f> 2 



+ v 2 p = (5) 



with those boundary conditions, so that orthogonality (which may easily be checked) and com- 
pleteness are guaranteed by Sturm-Liouville theory. We will later show that there is a one-to-one 
correspondence between these terms and the expectation value in irreducible representations of 
the gauge group; this is why we have picked this form. This means that any distribution which 
does not vanish at the boundary does not have an interpretation in terms of the underlying group 
theory. The large area limit is already visible as the a\ term. 

As a final piece of preparatory formalism, we also calculate the average Wilson loop in the 
fundamental representation from the Fourier series; this is 



(Wo) = 



p A (Wo)W dW = f p A (W ) sin t cos £ dcf>/2 

1 J-TT 1 1 



(1/4)/ (6 lS in^)sin^#= (6) 



since by orthogonality only the first odd term in the expansion survives. Later we will generalise 
this formula to all representations. 

3 The uncorrelated loop approximation 

We now claim that due to gauge invariance we can always write an individual Wilson loop (the 
qualification is important) chosen from the distribution as 



W = Wl + Jl - W$a.e r , (7) 
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Figure 2: Multiplication of adjacent Wilson loops. We can move the opening using a gauge 
transformation and hence create loops of any size in this fashion. 

where e r is a unit vector in a random direction in the 3-space inhabited by W. This is essentially 
a consequence of Elitzur's theorem. To see this, consider any Wilson loop of a given shape with 
its opening at a particular point on a lattice. If the distribution of W is not as shown, then there 
must be some preferred direction(s) for W. However, we can perform a local gauge transformation 
at the site of the opening. It is easy to see that this leaves Wo invariant, while rotating W to an 
arbitrary angle. As we have complete freedom to do this, with the single exception that rotations 
of Wilson loops with their openings on the same site are correlated, we can move the preferred 
direction of W. (We can avoid the exception here by simply considering only sets of translated 
Wilson loops without rotations, which gives a perfectly valid distribution.) Thus there is no 
preferred direction after all, and the distribution is the one claimed. 

We now take two Wilson loops, W 1 (^4) and W 2 (A), of the same shape and area A, but chosen 
in such a way that their opening lies on the same site, fig. ^. It will also be convenient to choose 
rectangular loops with a side in common traversed in such a way that this cancels out in the 
product of the loops. Thus we produce a Wilson loop of twice the area. The Wq element of this 
loop, by the standard SU(2) multiplication rule, is: 

W (2A) == W ( }(A)W 2 (A) - (1 - <(A) 2 ) 1/2 (1 - W 2 (A) 2 ) 1/2 cos£, (8) 

where cos£ = e^.e 2 .; £ is the angle between the two unit three vectors in the directions W 1 (A) 
andW 2 (A). 

The correlations between the diagonal elements Wg(A) in eqn. (||) together with the angle £ 
contain all the information about correlations between adjacent loops; without them we simply 
have loops from a random distribution multiplied together. If we were to go on and consider larger 
loops still, we would have to consider the cumulative effect of these correlations on the larger loops, 
and the expressions would rapidly become unmanageable. 

Instead, we make the approximation that both the angular and the diagonal correlations are 
negligible for loops larger than some area A — Ad- We shall later see explicitly that this gives the 
leading strong coupling value for (Wo), although in our case expanded in terms of smaller Wilson 
loops rather than the coupling itself. That such behaviour is expected is due to the random nature 
of the choice of loops from the distribution. 

A random distribution for cos£ is simply a flat distribution, which follows from the fact that 
the usual measure for the polar angle on a 2-sphere (here the one containing e* ) is sin £ d£ = d cos £. 
We are assuming, however, that the loop distribution p(Wo) itself has not reached the large area 
limit. The formula (j|) can now be used to multiply loops of different sizes, and we can build 
up larger Wilson loops by repeated application. Moving the opening in the Wilson loop simply 
corresponds to a gauge transformation, leaving the only remaining variable Wq invariant, so there 
is no difficulty in constructing Wilson loops of any area which is a multiple of Ad by this method. 

Our procedure is to take the measured Monte Carlo distribution of Wq on loops of area A = Ad 
as our input and hence use the approximation to calculate the distribution of Wilson loops for 
larger areas. In our examples, we shall in fact start from the plaquette, i.e. take Ad = 1 in 
lattice units and the measured plaquette distribution, but this is not forced on us. This method is 
obviously something of a hybrid, since the initial distribution includes all correlations, but it will 
allow us to perform explicit calculations using the loop distribution. Further, we will be able to 
predict the behaviour, starting from the appropriate expectation value for the plaquette, in every 
representation of the gauge group. 
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Calculation of the product distribution 

Suppose that the probability distribution of Wq for loops of the initial size Ad is given by 
PA d (Wo)dWo, normalised so that the integral between Wo = — 1 and 1 is unity. To produce 
the distribution of larger loops we must integrate over the smaller ones, using the equation (^) as 
a constraint. It will be helpful at this point to change to angular variables: take Wo(2A^) = cos 9, 
where < < ir, and likewise W^Aj) = 6i for < Oi < it. The new distribution is given by 

P2A d (W )dW = \ [ [ PA d ( C0S Oi) PA d ( C0S 62) sin 0^0! sin 6 2 d9 2 sin £d£„ (9) 



where the extra factor of 1/2 in eqn. (£)]) has appeared because in covering the full range of 6\ 
and 62 we cover the range exactly twice. This can be seen by simultaneously exchanging 
61 <-> 7r — 6\ , 82 <-> 7r — 62 in eqn. (|l0|), which leaves Wq invariant. 
The range in this notation is taken from eqn. (||) which becomes 

Wq = cos = cos 6*i cos 62 — sin Q\ sin O2 cos £. (10) 

Note that with the correct normalisation of ^^(cos^), the new distribution will also be cor- 
rectly normalised, i.e. / P2A d (Wo)dWo = 1, regardless of the form of p; this is nothing more than 
a basic consistency condition for the integrals we are performing. 

We remove dt; from eqn. (||) using eqn. ([To|), sin#i sin #2 sin£G?£ = dWo, giving 

P2A d (cosQ) = - I I pA d (cos6 1 )pA d (cos9 2 )d6 1 d92 (11) 

^ J J constraint 

with the constraint of eqn. ([l0]). The considerable simplification due to the approximation that 
cos£ has a random distribution is now evident. 

The explicit evaluation of the integral is rather uninteresting. We shall simply point out that 
it is easiest to change variables to the </> used in the Fourier series ([|), so that W = cos 6 = sin 5. 
We then insert the Fourier series for the initial distribution pA d into the integral, and after some 
even more tedious algebra best left to computers we obtain, 

P2A d (sin-) = K — A f^cos(n --)$ + -£ sian*) . (12) 

n— 1 ^2 ' 

Both odd and even sets of terms remain separate; as we have already hinted, this is no coincidence. 

The relationship between b\ and the loop expectation value was given in eqn. (6). We now read 
off the coefficient of sin $ which is 7rfe^/4, hence the new expectation value is (wbi/A) 2 . Because the 
terms in the expansion remain separate, we can immediately extend the result to N multiplications 
and a loop size A = N + 1, 

(W (A)) - (nh/Ay+ N - W (1) A , (13) 
an area law behaviour with string tension — log(Wo(l)). 



Eqn. (13) is equivalent to the lowest order strong coupling result, even though in this case no 
mention of the coupling has been made; indeed, it is renormalised due to the fact that we start 
from the physical plaquette rather than the strong coupling prediction, although the behaviour 
for larger loops in terms of the plaquette is the same. 

In fig. [?[ we compare distributions from uncorrelated loops with the measured values taken 
from fig. 0; as already stated, we have started from the measured plaquette distribution to generate 
the uncorrelated loop distributions. The multiplications in this case were done by Monte Carlo 
integration of the rule (^), which was rather more convenient than applying the analytical formulae, 
thus accounting for the noise in this distribution. 

Not surprisingly, we immediately see the measured shape start to diverge from the uncorrelated 
loop shape. So it is not at all clear this simple approximation is useful in explaining the physics. 
However, later we will see that the terms in eqn. (|l^) correspond to irreducible representations of 
SU(2) and that the individual terms separately obey an area law. Thus the calculation has the 
useful effect of demonstrating the leading strong coupling behaviour in all representations. 
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Figure 3: Comparison of Wilson loop distributions from Monte Carlo and those from the uncorre- 
latcd loop approximation (where the noise is clearly visible). The measured plaquette distribution 
has been used as a starting point for the latter. 



4 Centre dominance 

As explained in the introduction, the centre dominance picture applied directly to Wilson loops 
has been shown to be quantitatively successfully, indeed in very good agreement with the full 
values 0. In this approach any Wilson loop obtained from Monte Carlo results in SU(2) is 
approximated by its sign. It should be pointed out that in a weak sense it is not surprising that 
this reflects aspects of the dynamics: if one divides p(Wo) into odd and even parts, only the odd 
parts can contribute to (Wo) because that comes from the first moment of the distribution. We 
shall return to this point below. The result that the sign alone contains the essential dynamics is 
stronger and deserves examination. 

In terms of the distribution, the centre dominance picture requires that we count — 1 for a loop 
with — 7r < < and +1 for < <\> < ir with Wo = sin(0/2), i.e. a step distribution. We use the 
Fourier expansion of the Wilson loop distribution in eqn. (Q) and obtain for the expectation value 
in this picture, 

/w z(2) x f n (/>d<j) f° (f>d<p 

(W W )=/ PCOS-y- / PCOS — 







l^/f (f> f° <t> \ ^ 

— y I b n sin n</> cos — d(p — I b n sin n0 cos — -d<fi I = / 



(W) 



2 ^ Vo Y 2 - J_ n " 2 V ^ n2 - 1/4 

The series here is to be compared with the exact value nbi/4. The important difference is 
the presence of all the higher b n , with decreasing coefficients, in the centre dominance formula; 
without these the results would differ only by an overall factor of 37r/16, which vanishes in the ratios 
required for the heavy quark potential (see eqn. (p6|), below). Why the value in eqn. (14) gives the 
correct value can therefore by clarified by looking at the values of b n in the real distribution. Some 
support for this result is given by the uncorrelated loop picture, eqn. (|l2|), where in doubling the 
area of a loop the new coefficients b n are suppressed by a factor 1/n, which in combination with 
the factor n/(n 2 — 1/4) in eqn. (14) shows that higher terms become progressively less important 
both for large A and for large n. If the other b n are already smaller than b\ for small loops this 
effect is enhanced because of the powers of b n involved. 



7 



1 



0.1 
0.01 

c 
_Q 

0.001 



0.0001 



1e-05 

2 3 4 

Coefficient n 



Figure 4: The ratios of Fourier coefficients |& n /&i| plotted against n for the three smallest loop 
areas. The squares show the same values in the uncorrelated loop approximation. 

We have calculated the coefficients b n by binning Wilson loop data from Monte Carlo simula- 
tions and performing the integral for the Fourier coefficients, 

i r 

b n = — p sin n<bd<j) : (15) 

numerically on the data at the end of the run. Errors have been estimated by averaging over values 
of b n obtained in this way from several runs. The results here use 500 bins; we have compared 
this with the result from 200 bins in order to check that the discretisation error due to the finite 
number of bins is small. This error increases with the coefficient n as usual when sampling higher 
frequencies. 

In fig. [I] we show |6„/&i| for the three smallest areas (given in terms of the lattice spacing a) 
plotted against n on a logarithmic scale: I&4/&1I at area three is omitted because of large errors. 
The data is from a 12 4 lattice at /3 = 2.5. Plotted in this fashion, the data shows two things. 
First, b n drops off faster for larger n. Second, the rate at which this happens increases for larger 
loops, as witnessed by the increased slope. Even by area 3a 2 the value for n — 3 is down by a 
factor over a hundred on n = 2. 

Thus the 61 term is the only one contributing to the centre dominance result, eq. (14), for 
loops of even moderate size and hence the result is the same, up to a trivial normalisation factor, 
as that from the full Wilson loop for all medium and large loops. 

In fig. |], we have also shown the ratios for n = 2 and 3 using the uncorrelated loop approxi- 
mation, for comparison. As we have used the value in the plaquette distribution as the starting 
point for the higher coefficients, the points for area la 2 are chosen to be the same as in the full 
case. However, the trend thereafter is similar, too. It has come from two places: first, the intrinsic 
behaviour of the model seen above, where higher b n are more highly suppressed in larger loops; 
secondly, the fact that the initial coefficients b n for higher n are already smaller and so powers of 
higher b n disappear faster: only the first of the two is an explicit prediction of the model. 

Behaviour in higher representations 

It turns out we are able to explain the results in the previous section as the effect of higher 
representation loops. As an alternative to eqn. (14), where we expanded the expectation value in 
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the centre dominance picture in terms of Fourier coefficients, we could equally well have performed 
an expansion in terms of the characters of various representations; our conclusion would then have 
been that only the lowest half-odd- integer representation contributed for large loops. In fact, the 
expansions are the same up to constant factors. 

To show this, we shall define the character of the representations in terms of the corresponding 
traces of the Wilson loop as 

Xm(W ) = (m+l)Tr?W, (16) 

where the representation has isospin m/2, so that the Tr operator is defined to contain the normal- 
isation factor usually used in Monte Carlo calculations. The standard recurrence relation between 
the characters for SU(2) is 

Xm+l(Wo) = Xm(Wo) X l(Wo) - Xm-l(Wo). (17) 

We shall now show that using our variable (f>, defined in terms of the fundamental representation 
as Tr x /2 W = Wo ee sin ^, the characters for odd and even m can be written, 

(w\-t i^-i cos ( n ~i)^ 

X2{n-\){W( S ) = (-1) -g , 

cos | (18) 

X2n-l(W ) EE (-1)" T . 

cos ^ 

It is to be assumed that the value of the character at 4> = ±ir is found by taking the limit, although 
we have already noted that the distribution must vanish there from group-theoretic considerations. 
Note that Xq(Wq) = 1 and Xi(Wo) ee 2sin| as expected. Now we must prove that the following 
versions of eqn. ([17]) hold, 

X2n(W ) + X2n-2{W Q ) = X2n-1 (W )Xl (W ), 
X2n+l(W )+X2n-l(W )=X2n(W )xi(W ). 



This can easily be done by insertion of the formulae (|18|) into the respective left hand sides, whence 
standard trigonometric relations give the required results. Hence by induction the equations (|l8| ) 
are true for all n > 1. 

Now we work out the expectation values of the traces using the distribution. For the represen- 
tations with j = n — 1/2, n > 1, 

(Tr n _i W) = —J_ p(Wb)x2n-i(W )cos|-^ = 1 '— " , (20) 

after inserting the Fourier expansion of eqn. (|j) and the second of eqns. (|lq). Likewise for the 
representations j = n — 1, n > 1, we have, 

^-^ = W (21) 

where n = 1 and 2 correspond to the identity and adjoint representations respectively. Hence, as 
claimed, there is a one-to-one correspondence between the Fourier terms and the representations 
of the gauge group. (The basic point about the relationship between the sign of the loop and the 
character expansion was made in ref. p]].) 

This means that the behaviour of the coefficients seen in fig. [l| is entirely consistent with the 
Casimir scaling hypothesis, recently discussed in the context of centre dominance fl7|| , which notes 
that higher representations j = m+l/2 have a string tension roughly proportional to the quadratic 
Casimir operator, 

Kj w K c x j(j + 1). (22) 
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A larger string tension causes a stronger area law fall off, and hence from eqn. (|2(j) the Fourier 
coefficients b m decay faster in the loop area. Combined with eqn. (14), it can be seen that Casimir 
scaling implies centre dominance in the j = 1/2 Wilson loops. This result would be rigorous if 
within Casimir scaling one included the assumption that all loops, not just those in the region of 
linear confinement, were suppressed by a similar factor. As the ratio of string tensions between 
j = 3/2 and j = 1/2 is already 5 in eqn. (p2|), it is not surprising that centre dominance holds so 
well. 

Combining these results with the uncorrelated loop calculation in eqn. ([l2]) , we see that in the 
strongcoupling limit all half-odd-integer representations have an area law behaviour similar to 
eqn. @, 

(Tr n _,W(A))=(tl)l^y. (23 ) 
Indeed, we have also an area law in the integer representations, 

This is perhaps more of a surprise, as one often encounters the statement 'strong coupling predicts 
a perimeter law behaviour in the adjoint representation'. The point is simply that the random 
plaquette model we have used sees only the planar contribution to the adjoint Wilson loop Q The 
perimeter law term comes from a closed tube of plaquettes surrounding the loop, to which the 
planar plaquettes do not contribute in this approximation. In real loops one expects both to be 
present, hence at large R and T the perimeter law will eventually dominate. We shall later return 
to the question of centre dominance in adjoint loops. 



Centre projection and positive plaquettes 

In the last few years, detailed investigations of the importance of Z(2) degrees of freedom have been 
made by projecting the full SU(2) links to Z(2). One then makes plaquettes from these projected 
links, creating so-called projection vortices j|, a planar Wilson loop is simply a product of these 
over the minimal area. 

As we have shown (see for example eqn. (6)), any observable contributing to the Wilson loop 
expectation value must come from the odd part of the distribution, so that it changes sign when 
the relevant loop does. It is natural to suppose the projection vortices are related, by some effect 
of the projection, to signs of Wilson loops, since that is where the effect of the centre appears 
physically. If this is the case, then the observation that no vortices means no area law is trivial: 
as we saw in discussing the loop distributions, no odd-sign behaviour means, in analytical terms, 
no contribution to the real expectation value. This does not mean that the relationship between 
the vortices and confining behaviour is trivial — the claim is that the physical string tension can 
be seen at particularly short distances in centre projection — merely that the presence of both 
positive and negative signs is necessary to observe the dynamics. 

To make this a little clearer, consider the graphs showing the expectation value of Wilson loops 
separated according to whether they contain odd or even numbers of projection vortices in, for 
example, fig. 8 of ref. Q]: the Wilson loops separated in this way show no confining behaviour. 
Indeed, the fact that the Wilson loops go to zero at all is because the projection vortices do not 
exactly correspond to the signs of the full loops: suppose we look instead at the latter. As we have 
already seen, the distribution approaches (J3j) for large loops, and if we consider the part of the 
distribution with a positive value of the Wilson loop, we can see that large Wilson loops approach 
a constant value, 

(\W ol \) = J \Wo\pi(W )dW = ^ J | Wo 1(1 - W 2 ) 1/2 dW = (25) 

1 I am grateful to Jeff Greensite for making this point clear to me. 
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and the negative loops taken alone correspondingly approach the limit —4/ (3tt). This is a property 
of the even part of the distribution, unrelated to the odd part where (Wo) is to be found. (However, 
the adjoint string tension lives in the even part of the distribution: see below.) 

The relationship between thick vortices and the sign of the loop is clearer in the positive 
plaquette model Jl2| ]. Here, plaquettes are constrained to have a positive value: this can be 
implemented by a simple accept/reject test in Monte Carlo updating. This model appears to 
have all the right properties to be an alternative lattice regularisation of SU(2) [|l3|. There are 
no thin vortices, which are simply lines of negative plaquettes, so that through the commutativity 
property of the centre any Wilson loop with a positive (negative) sign must contain an even (odd) 
number of thick vortices. 

We have investigated the Wilson loop distribution in this theory, and show the results in fig. |5|. 
The coupling here is (3 = 1.9, found to be in the scaling region for this theory The approach 
to the large loop region is, unsurprisingly, very similar to that shown in standard SU(2) in fig. [l[ 

The sign of the Wilson loop here counts thick vortices exactly and so the expectation values of 
Wilson loops with only even or odd numbers of thick vortices must tend to ±4/37r, irrespective of 
whether the vortices are responsible for the dynamics, i.e. irrespective of centre dominance. Our 
point here is simply that one must consider loops with odd and even numbers of vortices together, 
as separately they can have no interpretation in terms of (Wo). 

We can also use the positive plaquette model to examine the question of whether the results 
from the Z(2) part agree with the full results when only thick vortices are present, or whether, on 
the other hand, both thick and thin vortices are required to produce the full result. Fig. |^ shows 
the potential between a heavy quark and an antiquark calculated from 150 configurations on a 
16 4 lattice at (3 — 1.8 for both the full results and those using only the signs of the Wilson loops. 
This figure corresponds to figures 7 and 8 of ref. ]t| which show results in standard pure SU(2); 
our method is the same, i.e. V(R) is calculated from 

V(R) = lim - log W ( R > T+1 ) (2 6 ) 

with the signal improved by using fuzzing on spatial links only — this does not affect the presence 
of vortices in the (R, T) plane. To enable a comparison with the standard case we have also made 
a fit to the form 



V{R) = V Q + KR-e/R] (27) 

no correction for lattice artifacts has been made and the single correlated fit included all data 
with T > 3 and R > 2.5 for the full loops only (i.e. not using the Z(2) projected values). We 
obtained V = 0.572(6), K = 0.0410(9) and e = 0.261(11), corresponding to slightly below 2.5 
for standard SU(2). 

It is clear from fig. ||that, as with standard SU(2) in ref. 0, the Z(2) part carries the physics; 
a similar picture (not shown) holds for (3 = 1.9 where the curvature is more pronounced. Thus 
centre dominance is present across the whole size range, from the Coulomb to the confining region, 
works even with no thin vortices present. This is something of a relief, as the division into thin 
and thick vortices is an artifact of the cut-off. To see this, consider taking a negative plaquette, 
i.e. part of a thin vortex, and halving the lattice spacing: on the new lattice, it is no longer clear 
whether the loop corresponding to the former plaquette is (say) part of a thick vortex surrounding 
four positive plaquettes on the new lattice, or is negative by virtue of one of the those plaquettes 
being negative. In other words, the new lattice does not preserve the division into thin and thick 
vortices of the old one. 



Confinement from non-interacting vortices 

We will show that a random distribution of non-interacting thick centre vortices, with no other 
dynamics, can give an area law in SU(2) with a string tension of about the right magnitude. This 
is the simplest possible treatment, although it is very much in the spirit of Nielsen-Olesen flux 
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Figure 6: The heavy quark-antiquark potential from the positive plaquette model at (3 = 1.8, in 
lattice units. The dashed line is a three-parameter fit. 
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Figure 7: The random non-interacting vortex model, as described in the text. 



vortices [l5| ; random vortices and their confining nature are not new and were considered for 
example in |16| . However it is now possible to make a quantitative test of the idea using results 
from centre vortices. A similar (but not identical) calculation appeared in ref. fl7j| ; the version 
here has the advantage that it makes no reference to the lattice, even if the underlying assumptions 
are too simplistic for it to be a true continuum theory. 

We assume that the value of a Wilson loop of area A depends only on the number of vortices 
inside, i.e. Wo(^4) = Wc{— 1)™, where n vortices pass through the loop Wq and Wc depends only 
on the coupling: this is just the centre dominance proposal, which we have seen works well in 
practice. Although the assumption in this form does not fit in with the Wilson loop distributions 
we have shown, the proposal is that the sign alone determines the dynamics, so that the effect of 
the rest of the distribution on (Wo) can be absorbed into Wc- 

We also assume that the vortices are independent from one another, so that any vortex can 
pass at a random point through a given area regardless of the presence of another: this second 
assumption is clearly an oversimplification, as there is indeed interaction between the vortices [ fl9| . 
We are also assuming the vortices are infinitely thin, so that there are no edge effects; this, too, is 
unphysical, as we have already stressed that the vortices must be thick to survive in the continuum 
limit. 

Finally, we suppose that the only relevant physical quantity is the mean density of vortices 
passing through unit area which we denote by pa- In addition, we are using the fact of translational 
invariance so that there are equal probabilities for a vortex to pierce any regions with the same 
area. We are also essentially ignoring the fact that the vortices appear from the gauge dynamics. 

Consider an area very much larger than A, A' = AN for N ^S> 1, which encloses A and through 
which exactly paA' vortices are assumed to pass, randomly distributed across the area: a three- 
dimensional slice is shown in fig. [?|. For each vortex independently, there is a probability 1/N that 
it lies inside A. The probability for n of them to lie in the area A is then given by a binomial 
distribution, 



PA-A'—n 

Pr(n vortices in A) = ) ^- j - - ) . (28) 

The contribution to the Wilson loop is 

PaA' , / i \ n / 1 \PAA'—n / r.\p A AN 

(W„(A)) = W C g ( W „ )("!)" (i) (l-i) -*,(.-£) . 
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Letting iV — > oo, and noting that the limit as N tends to infinity of (1 — 2/N) N is e 2 , we have 



which is an area law with string tension K = 2pA- This result was obtained in a different 
approximation by ref. ||; in this case we stress that its validity is based only on the notion of 
randomly distributed, non-interacting and (unfortunately) thin vortices. 

A string tension of K — (440MeV) 2 would hence require pa ~ 2.5 fm~ 2 . Determining this 
quantity in a gauge invariant fashion is difficult, so we shall use the value calculated directly from 
Z(2)-projected vortices in ref. |L8| : pa — (1.9 ± 0.2) fm -2 , close to what we require. Further, 
the same value of the string tension was assumed in that calculation, so in fact the ratios of the 
two numbers are independent of the experimental value of K. Hence the random vortex model 
predicts a value of K about three quarters of the measured value (giving for y/~K some 85% of the 
measured value), the difference corresponding to three standard deviations. This is surprisingly 
good considering the model in question contained very little mathematics and did not even take 
into account the extended nature of the vortices. That the centre vortices are physical is supported 
by the result of ref. Jl8]] that the distribution is renormalisation group invariant. 

The adjoint representation 

It has long been argued, and has to some extent entered the folklore, that the observation of a 
string tension in the adjoint representation of SU(2) cannot easily be explained by centre vortices. 
This claim is based on the point that the centre Z(2) of SU(2) is not seen by adjoint fields. 
This is, of course, correct, but this Z(2) is associated with thin centre vortices, which are not 
expected to be involved with confinement in the continuum limit — and as we saw from the 
positive plaquette model are not necessary even for lattice physics. The surviving symmetry of 
the adjoint representation is SU(2)/Z(2), which is just that required for thick centre vortices. The 
mathematics of confinement relates to the fact that with SU(2)/Z(2) one has two types of path in 
the gauge manifold, one contractable and one not; the latter are closed due to the identification of 
opposite ends of a diameter of the manifold. This distinction survives in the adjoint representation. 

In the adjoint representation, up to constants irrelevant to the argument, the action depends 
on the square of the plaquette value, and the negative trace of the loops is not seen. However, 
the fact that one cannot measure the signs of the loops in no way blurs the distinction between 
the topologically different windings of the field. One has simply lost a useful but — if the theory 
of thick centre vortices is correct — physically unnecessary label, namely the association between 
the sign of the loop and the presence of a vortex. 

This last remark explains the difficulty with projection vortices, where factors of —1 presum- 
ably correspond to a thick centre vortex with an increasing admixture of thin vortices for stronger 
coupling. Such vortices certainly are not present in the adjoint representation, making investi- 
gation of the mechanism considerably more difficult, but the underlying topology of the gauge 
manifold which forms the basis for confinement in the theory is still present. 

The point made here is slightly different from the suggestion |17| for reconciling the adjoint 
string tension with the vortex picture involving sign flips. However, the finite width of the vortices 
enters in both cases: in that mechanism it was the origin of string behaviour in higher representa- 
tions, while here we have stressed the SU(2)/Z(2) nature of the loops, which must be spread over 
a large area so as to be present in the continuum limit. Thus there may be a connection between 
the two. 

We can show the correlations between the behaviour of adjoint Wilson loops and the presence of 
vortices directly. We have run a simulation on our 12 4 lattice at (3 = 2.5 in which the adjoint loops 
are separated according to whether the fundamental loop was positive or negative, i.e. whether 
there is an even or odd number of vortices passing through the loop. This is just the process we 
criticised above for the fundamental case; however, here the sign of Wo itself is not seen by the 
adjoint trace, 



(W (A)) =W c e 



-2p A A 



(30) 



Wr 



radj 



(4W 2 - l)/3, 



(31) 
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Figure 8: Logarithms of adjoint Wilson loops plotted against area. Both the complete data, and 
those loops associated with a positive fundamental loop, are shown. 



so that any correlation with the sign involves the SU(iV)/Z(JV) topology associated with the thick 
vortices. Indeed, the adjoint loop probes the even part of the distribution p(Wo) whereas, as we 
have made clear, the expectation values of fundamental loops and their signs only involve the odd 
part. 

In fig. ||, we show the logarithms of the adjoint Wilson loops, both for the complete data 
and those with positive fundamental trace, once more at (3 = 2.5. There is clear sign of an 
area law behaviour in the former, i.e. an adjoint string tension. The latter have a consistently 
shallower gradient, in other words the string tension is significantly lower in those loops with only 
even numbers of vortices. A perfect or near-perfect correlation is not to be expected since from 
eqn. ( pi] ) that would require all loops with Wo < also to have \Wo\ < 1/4, and the other way 
around for positive Wo- 

Phcnomenologically, the source of the correlation is not hard to see from the fundamental 
distribution in fig. bearing in mind that the expectation value we want is just eqn. ( Jjjj ) integrated 
over this distribution. For smaller loops, the distribution peaks at positive Wo and tails off towards 
negative values. Therefore, large absolute values of Wq, and hence positive contributions to W J , 
tend to come predominantly from positive fundamental loops. However, the source of the adjoint 
string tension itself is rather less easy to fathom by this method, and for that we appeal to the 
strong coupling result of eqn. ( pi]) with n = 2. 

As this work was being completed, ref. |20|| appeared, showing similar correlations between the 
presence of vortices and the adjoint string tension. 



5 Summary 

The major results of this paper are as follows: 

• We considered the distribution p(Wo)dWo of Wilson loops in SU(2); for loops larger than 
the correlation length of the field p(Wo) reaches the limiting form (2/7r)^/l — Wq. 

• We applied the Fourier decomposition in eqn. (Q); the expectation value of the loop is 
(Wo) = tt&i/4, in other words it only probes the first odd coefficient of the distribution in 
this parametrisation; the even portions of p(Wo) do not contribute. 
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• In the simple approximation where all Wilson loops are uncorrelated, and we use the mea- 
sured plaquette distribution as input, the coefficient b\ and hence (Wo) obey the strong 
coupling behaviour, eqn. (|l3|). Larger coefficients b n are suppressed by an additional factor 
1/n for each multiplication by a Wilson loop of the same size. This approximation also 
shows an area law behaviour for the adjoint representation. 

• The Z(2) expectation value contains all b n , not just b\. We found that for n > 2 the 
coefficients were strongly suppressed in Monte Carlo results, showing the validity of centre 
dominance. 

• The Fourier terms correspond, term by term, with the expectation value of the Wilson loop 
in the gauge group representations. The faster decay of higher representation loops therefore 
gives an explanation for centre dominance. 

• We stressed that care needs to be taken when looking at results in Z(2)-projection (whether 
the links are projected or we simply look at the sign of the loops): the expectation values 
of the negative and positive parts separately are trivially seen to be non-confining, from the 
basic properties of the distribution. 

• We showed by considering the positive plaquette model that the absence of thin vortices, 
which are chains of negative plaquettes, does not change the centre dominance behaviour. 

• We showed that even a simple physical model of a random gas of centre vortices with no 
interactions gave three quarters of the observed string tension from the measured vortex 
density. 

• We pointed out that there is no essential difficulty in having thick centre vortices explain an 
adjoint string tension, as both the thick vortices and the adjoint representation are associated 
with an SU(iV) /Z(iV) symmetry — the theory of centre dominance does not directly involve 
the centre Z(7V), as has sometimes been incorrectly assumed. 

• We showed that there is a correlation between the adjoint string tension and the sign of 
the corresponding fundamental Wilson loop, even though the latter does not directly affect 
the adjoint loop: this makes it plausible that thick vortices are involved here, too. This 
correlation is easily seen by looking at the Wilson loop distributions. 

These results seem to show that there is indeed physics contained in the centre dominance 
picture. Of course, the results for other pictures of confinement need to be borne in mind. For 
example, a lot of physics has also emerged from the maximally Abelian monopole condensation 
(Abelian dominance, or dual superconductor vacuum) picture |plf ; see e.g. the proceedings of 
Lattice '97 p3| for recent progress. A way in which Abelian monopolcs might be related to centre 
vortices has been proposed || . It is interesting to note that both models are based on topological 
mechanisms using homotopy classes, those of SU(2)/Z(2) or U(l), not manifest in the full gauge 
group, and one can speculate that there is a deeper connection. There are of course other models, 
too, such as the anti-ferromagnetic vacuum p3| . The arguments for the fundamental importance 
of Z(N) in SU(N) given by 't Hooft twenty years ago ]2ij are still valid and should be born in 
mind, but the debate looks set to continue. 

For further progress, one needs a better way of asking which properties are the most funda- 
mental. Work is currently in progress by another group to try to understand the correlations 
between the monopole picture and instanton effects |25| ; centre dominance is presumably another 
ingredient which needs to be taken into account. How far our results here really reflect dynamics 
involving the centre of the gauge group can presumably be clarified by looking at SU(3). Recent 
results [^6| suggest the corresponding behaviour of Wilson loops for the centre Z(3) is indeed found 
in the SU(3) case. 
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